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In doped two-leg spin ladder systems, holes are expected to form charged bosonic pairs. We study 
charge transport in this system in the temperature range where the bosons can be described as weakly 
interacting quasi-particles. We consider boson-phonon and boson-impurity scattering processes. We 
suggest that due to the loffe-Regel resistivity saturation mechanism, the resistivity may exhibit a 
local maximum at intermediate temperatures. We propose that this may explain a similar feature 
found in recent experimental results. 
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I. INTRODUCTION 

One of the most exciting advances in the field of high 
Tc superconductivity and in the theory of strongly cor- 
related electrons in the last few years has been the ex- 
perimental study and the theoretical understanding of 
pure and doped quantum spin ladders |0| . In contrast to 
the one dimensional (ID) antiferromagnetic (AF) Heisen- 
berg chain, where the ground-state has an infinite corre- 
lation length, and exhibits a slow decay of the spin cor- 
relations, even -leg (undoped) ladders have spin-liquid 
ground-states with purely short-range spin correlations. 
The exponential decay of the spin-spin correlation is pro- 
duced by a finite spin gap, i. e., a finite energy gap to the 
lowest S' = 1 excitation in the infinite ladder. Even-leg 
ladders are the realization of the Resonant Valence Bond 
(RVB) quantum liquid of singlets proposed by Anderson 
in the context of 2D S — \/2 AF Heisenberg systems 

Two- leg S = 1/2 ladders are found in vanadyl 
pyrophosphate {VO)2P207 and in some cuprates like 
SrCu203 Spin susceptibility, neutron scattering, 

NMR and muon spin resonance measurements are con- 
sistent with a spin gap. 

Even-leg ladders are especially interesting under dop- 
ing because theory has revealed that holes should pair to 
form bosons, in a relative "d-wave" state, with a super- 
conducting ground-state at finite doping Copper- 
oxide compounds exhibiting one dimensional features 
with both isolated Cu02 chains and CU2O3 ladders 
(pairs of Cu02 chains linked by oxygen atoms between 
the coppers) have recently been found to be supercon- 
ducting d,^. In the compound Sr2C ai2Cu2iO n , a cop- 
per valency of~ 2.2, i. e., a hole density of 0.2 per 
ladder-Cu is estimated 1^. NMR and conductivity mea- 
surements have been conducted on the latter compound, 
and the existence of the spin-gap is clearly established 
for pressures smaller than 29 kbar 

Thus, in the normal state of two-leg compounds the 
charge carriers are likely to be bosons, restricted to move 
preferentially in one dimension (ID). The proposal that 
CuO-planes in high-Tc cuprates are Bose conductors in 



the normal state has been extensively discussed in con- 
nection with Anderson's RVB model Q. However, it is 
still controversial. Besides, bosons in this context are 
strongly coupled to fermionic spinous via a gauge field 
0. The situation is clearer and simpler in the doped 
spin ladders, where the bosonic carriers are not coupled 
to low energy spin excitations. Materials in which these 
ladders occur are good candidates for true bosonic con- 
ductors. 

At low enough energy, a one dimensional bose gas can 
be described in terms of density fluctuations with an en- 
ergy spectrum which is linear in wave vector |p^ . A 
one dimensional bose gas can thus be mapped on a one 
dimensional Luttinger liquid of fermions, the transport 
properties of which are known [ pd| : in the presence of a 
random impurity potential, Anderson localization at low 
temperatures causes the resistivity to increase to infin- 
ity as a power law of the inverse temperature. At larger 
temperature a metallic behavior may step in, but no in- 
termediate maximum is expected. 

There is a characteristic energy scale above which the 
collective density oscillation picture ceases to be valid, 
and the bose gas can be described by the semiclas- 
sical quasi-particle picture ||T^. Assume that the in- 
teraction potential between bosons can be modeled by 
c5{xi — Xj), where c is a constant. Then the excita- 
tion spectrum remains that of independent quasi par- 
ticles for Ti^k"^ /2rac ^ cub, where nic is the boson mass 
and Tib is the linear boson density. Thus for ksT > crib 
the ID Bose particle transport phenomena can be de- 
scribed within a semi-classical approach. At low density 
and small interaction strength c, this temperature can be 
small. Below fcsT < crib the collective mode description 
of the bose gas should hold and the transport should then 
be described accordingly pd| ]. 

One of us pointed out a long time ago § that in Bose 
conductors the condition for the validity of semiclassi- 
cal transport theory - that the carriers mean free path 
is greater than their De Broglie thermal wavelength - 
is much weaker than that for ordinary (Fermi) conduc- 
tors, for which the mean free path has only to be greater 
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than the mean distance between carriers. It was sug- 
ested in Ref. j|] , on the basis of the loffe-Regel criterion 
that when the mean free path is smaller than De 
Broglie's thermal wavelength the resistivity saturates at 
a temperature-dependent value, which may be estimated 
by replacing the mean free path by De Broglie's thermal 
wavelength in Drude's formula for the resistivity. 

The purpose of this note is to examine the conse- 
quences of the semiclassical picture on the resistivity 
of doped quantum spin ladder systems, in the spin 
gap regime. The hole pairs are treated as structure- 
less bosonic particle: no account is taken of the d-wave 
like symmetry of their wave function 0|. In the fol- 
lowing we examine the temperature dependent resistiv- 
ity of a one dimensional boson liquid interacting weakly 
with phonons and with impurities. We propose that 
this straightforward analysis leads to a qualitative un- 
derstanding of the experimental behavior of the electri- 
cal resistivity with temperature in doped quantum spin 
ladders ||. 



II. RESISTIVITY 

We assume that the free boson states are described by 
an energy band e(k) the bandwidth of which for k per- 
pendicular to the ladders direction (c-direction) is small 
compared to ksT. For motion along the ladders we as- 
sume that the bandwidth Wc is large compared to fcsT. 
Under these assumptions we may approximate e(k) by 
e(k) = h^k^/2mc, where — tt/cc < kc < n/ac is the com- 
ponent of k along the ladders direction, ac is the period 
within the ladder and rric is the effective mass. 

The calculation of the resistivity follows the traditional 
Boltzmann equation approach, as described by Ziman 
. We consider the contributions from boson- impurity 
and boson-phonon scattering to the resistivity. 



A. BOSON-IMPURITY SCATTERING 

The expression for the impurity scattering contribution 
to the resistivity, pi is ]l3| , |l4t . 
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In Eq. (2.1) , u is the bosons drift velocity along the lad- 
ders direction, Vk = 9e(k)/(9k and is the transition 
probability, given by 



<5[6(k)-e(k')]4(l 
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where nj is the number of impurities by cm? and Vb-i is 
the boson- impurity interaction matrix element, assumed 
constant. 
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is the boson occupation number, and /i is the chemical 
potential. 

The relation between /i and the number of bosons per 
cm in the ladders, Ub, is obtained in the usual way, as- 
suming non-interacting bosons p5| . We find 



nhXc = gi/2{z) , 
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where Ac = jmckBT is De Broglie's thermal wave- 
length for motion along the ladders, z — e^^l^^T 



Carrying out the integrations in Eq. (2.1) we find that 
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where n is the number of bosons per cm^ and rig = 
z/(l — z) is the zero-momentum state occupation number 
P6| . Note that the T-dependence of pi come solely from 
the ng term. The latter is a monotonously decreasing 
function of T. The behavior of the impurity contribution 
to the resistivity with T is unusual because it increases 
as T decreases. This is a consequence of the bosonic ID 
character of the problem. 



B. BOSON-PHONON SCATTERING 

The boson-phonon scattering contribution to the resis- 
tivity, pp is given by [|3|, 

Jd^kd^k'd^Q[ik-k')-u)rvt'.Q 
PP = — ■ a;^^— ' (2.6) 
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where q is the transition probability. Assuming a de- 
formation potential-like interaction between the bosons 
and the phonons and a Debye phonon spectrum, with 
sound velocity s, this quantity is given by 

K[q = f «Q'5(k - k' - Q)^(e(k) - e(k') - nsQ) 

xriQuiil + ni,) , (2.7) 

where a is a constant that characterises the boson- 
phonon interaction and nq is the phonon occupation 
number. 

We write pp as 
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where kjj is Debye's wavevector, Qo ~ hskp)/kB is De- 
bye's temperature, 

Jo Jo 

(2.9) 
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n{x) = (e^ - and x±{v) = ^1^3^/21^ ± sQ/2 + 

g^r/VSmc - ^i|kBT. 

At low temperatures, T, Y+{rj) ~ Q^ry^/Smc 

and it follows from Eqs. (2^) and (2^) that 
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At high temperatures, T ^ and assuming that 
rricS^T lkB&£, ^ 1 and UcXc <C 1 we find that 
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C. RESISTIVITY SATURATION 



The expressions for p used above, Eqs. ( |2.l| ) and ( |2.6| ) 
, are valid as long as the mean free path for motion along 
the ladders £c is greater than Ac- For < Ac resistiv- 
ity saturation is believed to occur. An estimate of the 
saturation value of the resistivity psAT is obtained by 
substituting (.^ by Ac in Drude's formula p — rric/nq^T. 

Defining £c as — ^tt, whe re vt is the bosons ther- 
mal velocity, defined as vt = •\/2(e)/mc, where (e) is the 
bosons m ean thermal velocity. Using the boson distribu- 
tion Eq. (|2.3|) we find that 
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where 53/2 = J2T=i 
Thus 
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III. DISCUSSION 
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According to Eq. (2.13) psAT increases with increasing 
T, whereas pi increases with decreasing T. Thus, if the 
boson-impurity interaction is sufficiently strong, there is 
a temperature below which pi > psAT- Since, accord- 
ing to the loffe-Regel criterion the resistivity cannot be 
greater than psAT, saturation occurs and the resistivity 
curve follows psAT- Thus there is a resistivity maximum 
at a temperature To that can be estimated as that for 
which pj ~ psAT- This is the main result of this paper. A 
typical resistivity curve predicted by our model is shown 
in Fig. |l|. This figure also shows experimental resistivity 
data for SrxCai-xCu2404i at l.SGPa obtained in Ref. 

In Fig. d the contributions to the resistivity from 
phonons, impurities and from the saturation mechanism 
are normahzed as pi{T = lOOK) = psat{T = lOOA') = 
22mn ■ cm and pp{T = 200A:) = 5mfi • cm. These val- 
ues for pi and psAT are chosen so that the theoretical 
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prediction coincides with the experimental one. The pp 
curve is a smooth interpolation between the T <^ 9p) and 
T 9d behaviors, with 6d = 20QK. The value of pp at 
T — 60 = 200 AT is chosen for illustrative purposes only. 

We propose that the local resistivity maximum ob- 
served in the SrxCai-xCu2A04^i data at T ~ lOOA' and 
several pressures Q arises from the same mechanism as 
that for the ID boson model discussed above. Thus, 
we interpret the decreasing resistivity curves below and 
above the temperature where it is maximum as result- 
ing from resistivity saturation and from boson-impurity 
scattering, respectively. In order for this interpretation 
to be consistent it is necessary that, at T lOOA T, 
Pi ~ PSAT ~ 2QmVl ■ cm. 
this requires that mc/men ^ 
is the mass of the electron. 

We show in Fig. ^ the mean free path along the lad- 
der direction obtained from the SrxCai^xCu2AOiii data 
at 1.5GPa |^ using the Drude formula f or th e resistiv- 
ity and the definition of Ic given in Sec. II C . In Fig. 
^ the values of and n are chosen so that £c = -^c 
at T = lOOAT. This requires, as discussed above, that 
mc/meTi = 4.8 X 10~"'^^cm'^. We choose mc/m^ = 180 and 
n = 3.7 X lO^^cm"^, which is equivalent to 10~^ holes 
per ladder Cu atom. It is clear from Fig. ^ that there is a 
crossover around T ~ 100 AT from a regime of rapid varia- 
tion of £c with T, that according to our interpretation re- 
sults from boson-impurity scattering, to a regime of slow 
variation of 4 with T, from T - lOOA' to T - 50A'. The 
latter is consistent with the ideia of resistivity saturation, 
as discussed in Ref. jsj. We also note that, according to 
our proposal, the boson-impurity scattering contribution 
to £c can be estimated for T < lOOAT by extrapolating 
the high-temperature £c x T data to lower T values. As 
shown in Fig. ||, this leads to a rapidly decreasing £c xT 
curve. This further indicates that, within the semiclas- 
sical picture, resistivity saturation is bound to occur in 
this compound, independent of the precise value of 4/Ac 
at which it sets in. 

The numerical values for TOc and n obtained above are 
rough estimates. Neverthless, our picture requires values 
of mc which are too large compared to those obtained 
from theoretical estimates for pure two-leg ladders |l7| ] 
and values of n that are smaller than those reported in 
experiments (~ 0.2 holes per ladder Cu). However, the 
possibility of large effective mass enhancement due to 
boson-phonon interaction effects, neglected in theoretical 
considerations, cannot be excluded. Besides, the loffe- 
Regel criteria gives only very rough order of magnitude 
estimates for the values of £c/Xc and of p at saturation, 
so that the numerical values for and n used in our 
calculations may be uncertain by a large factor. 

For the semiclassical picture to be valid a second con- 
dition must be met: the cross over temperature 9 to 
collective bosonic behavior must be smaller than the 
temperature at the resistivity maximum. Our order of 
magnitude estimate for this relies on the fact that the 
hole pair spatial extension is S, so that the change in 
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binding energy AiJ of a pair when colliding with an- 
other one is AE ^ Eac/6, because two hole pairs can- 
not occupy the same rung. Since E ^ 200-?^, we find 
= cub/kB ~ Ent5{aJ6) - 20K. This value of 9 fur- 
ther supports our proposal. 

To conclude, we find that the simple minded semiclas- 
sical transport theory proposed here for the bosonic hole 
pairs in the quantum spin ladder system may account for 
the observed dependence of the resistivity on the temper- 
ature observed for SrxCai-xCu2404i. Our predictions 
can, in principle, be tested by changing the impurity con- 
centration, n/. According to our scmiclassical model, 
increasing nj should shift the resistivity maximum to 
higher temperatures and increase its intensity, whereas 
decreasing nj should shift the maximum to lower tem- 
peratures and decrease its intensity. 
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tions stated in Sec. Q , the semiclassical approach is not 
valid below this energy so that the fcc-integrals do not 
extend to kc = 0, but have a low fcc cutoff. We check 
that this condition is met in the numerical calculations 
performed here. 
[17] Steven R. White and D. J. Scalapino Phys. Rev. B 55, 
(1997) 6504. 

FIG. 1. Thick line: calculated resistivity versus tempera- 
ture curve for rric = ISOme and 10~ holes per ladder Cu. 
Squares: 5r2Cai2Cii2404i data. Dashed lines: individual 
contributions from phonons, impurities and the saturation 
mechanism (see text). 

FIG. 2. Mean free path versus temperature for 
Sr2Cai2Cu2404i at 1.5GPa . Dashed curve: extrapolation 
form high temperature data. 
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